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We have developed a novel system consisting of a superconducting single-electron transistor (S- 
SET) coupled to a two-dimensional electron gas (2DEG), for which the dissipation can be tuned 
in the immediate vicinity of the S-SET. To analyze our results, we have developed a model of 
the environment for S-SET/2DEG systems that includes electromagnetic fluctuations coupled both 
through the S-SET leads and capacitively to the S-SET central island. We analyze this model, 
treating the leads as finite transmission lines, to find the probability function P(E) for exchanging 
energy E with the environment. We also allow for the possibility of low- frequency fluctuations of 
the S-SET offset charge. We compare our calculations with measurements of SET conductance 
versus 2DEG conductance and find good agreement for temperatures > 100 mK, while unexplained 
discrepancies emerge for lower temperatures. By including the effects of charge averaging we are 
also able to predict the shape and evolution of I-V curves as the 2DEG in the vicinity of the S-SET 
is changed. 

PACS numbers: 74.50.+r,73.23.Hk,74.40.+k 



I. INTRODUCTION 

The effects of the electromagnetic environment on elec- 
tric transport has been a subject of extensive theoreti- 
cal and experimental interest in recent years. The rea- 
sons for interest are varied, as are the systems for which 
studies of the effects of the environment have been-gerp. 
formed. Recent interest in quantum computationtmocl 
has prompted interest in the effects of dissioation on 
decoherence rates in superconducting qubitsBtl Double 
quantum dots have been used to study the. effects of 
the environment on inelastic tunneling rates,Q and have 
been proposed as detectors o£ high-frequency noise pro- 
duced by mesoscopic .devices. □ Finally, interest in quan- 
tum phase transitionsO has prompted study of the effects 
of dissipation on superconducting systems such as thin 
filmst3'Ell and Josephson junction arrays.ll3 

Gjpjn the possibility of using a single Cooper pair 
boxclfca as a qubit and its similarity to the superconduct- 
ing single-electron transistor (S-SET) ,Ej it seems logical 
to use the S-SET as a model system for studying the 
effects of environmental dissipation on coherence and 
transport in small tunnel junction systems. Recently, 
there have been several such attempts, motivated by ex- 
periments in which a Josephson junction array was fab- 
ricated in close proximity to a two-dimensional electron, 
gas (2DEG) in an GaAs/Al^Gai-^As heterostructurejt 2 ] 
which can be used as a tunable source of dissipation. 
(While the effects of a mechanically tunable environment 
were earlier_studied in the macroscopic quantum tunnel- 
ing regimejij use of a 2DEG allows more flexible tuning 
of the environment over a larger impedance range.) In 
place of an array, we and the Berkeley group have instead 
used simila r fab rication techniques to couple an S-SET to 
a 2DEGJl3'tZlt3 In one instance, the focus was on trans- 
port at higher bias^i, in the regime of the Josephson- 
quasiparticle cycle.E-3 In the others, the focus instead was 



on th e low bias regime and the tunneling rates of Cooper 
pairs.l 16 H 18 l Theoretical work aimed at an explanation of 
the results of, the Berkeley group was undertaken by Wil- 
helm, et alS3 

The primary experimental difference between our own 
work and that of the Berkeley group lies in the way in 
which the environment is varied. The Berkeley group 
followed an approach developed earlier for the study of 
junction arrays,E£l using a gate on the back side of the 
substrate to vary the sheet density n s of the 2DEG, and 
therefore its resistance per square R sq . Such a change 
is global, and affects the 2DEG not only immediately 
beneath the S-SET, but beneath the macroscopic leads 
used to measure it as well. In our own work, by placing 
Au gates on the surface of the sample near the S-SET 
itself, we were able to vary the dissipation in the 2DEG 
locally, while. .Leaving 2DEG beneath the leads virtually 
unchanged .E§E] 

The Berkeley group compared their results to the the- 
ory of Wilhelm, et al., which predicted that within linear 
response the conductance Gset of the S-SET would scale 
with the ground plane conductance G*2D = 1/Rsq an d 
temperature T as G^/T". While the Berkeley group 
did observe power law behavior, their measured expo- 
nents were not in quantitative agreement with theory. 
Furthermore, the measured (3 depended on T and a on 
G2D, calling the scaling form into question. 

In our own work,ta we examined a somewhat more 
complex model for the environment than that considered 
by Wilhelm, et al. Specifically, we allowed for coupling of 
electromagnetic fluctuations to tunneling Cooper pairs 
due simultaneously both to the S-SET leads and to the 
2DEG in the immediate vicinity of the S-SET, which is 
coupled to the S-SET central island via a capacitance 
C*2D- By also allowing for averaging of the S-SET offset 
charge, we were able to obtain good agreement between 
our measurements and calculations. Here we examine 
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our model of the environment in somewhat more detail, 
provide additional experimental data which supports our 
earlier analysis, and also give some additional details of 
the calculation. 



II. SAMPLE AND ENVIRONMENTAL 
CALCULATIONS 

A. Sample Design and the Environmental Model 

1. Sample Design 

Our samples consist of an Al/AlO^-based S-SET fab- 
ricated in close proximity to a ,2I)EG formed in a 
GaAs/ALjGai-^As heterostructureE^I as shown in Fig. [j] 
below. We begin by fabricating six Au gates which can 
be used to deplete the electrons beneath them by appli- 
cation of a negative gate voltage V g . At the center of the 
Au gates we then fabricate our S-SET, as can be seen in 
the electron micrograph (expanded view in Fig. |l|). Note 
that as shown in the larger diagram the S-SET leads ex- 
tend over the 2DEG to macroscopic contact pads. For 
the vast majority of their length they are well away from 
the Au gates and the 2DEG beneath them is indepen- 
dent of V g . When we apply a gate voltage V g to all six 
Au gates, the electrons immediately beneath them are 
depleted, leaving a small pool of electrons beneath the 
S-SET. This pool is connected to the rest of the 2DEG 
(held at ground) only by two quantum point contacts 
(QPCs) with conductances l/i?Qpc (assumed equal) as 
shown in the micrograph. It is also capacitively coupled 
to the S-SET island through a capacitance C2D as shown 
in the lower right inset to Fig. [l| When all six Au gates 
are energized as described above, we say that the elec- 
trons are confined in the "pool" geometry. We do not 
refer to the pool as a quantum dot since for these exper- 
iments the QPCs are sufficiently open that no Coulomb 
oscillations are detected in the pool and discrete energy 
levels have not formed. 

Because the Au gates can be biased independently, we 
can also apply V g to only the four outermost gates which 
form the QPCs. As before, the electrons beneath the S- 
SET are coupled to ground through the QPCs. In addi- 
tion, however, they are now coupled through a resistance 
Rstr to two large reservoirs of electrons located between 
the four outermost gates, as can be seen in Fig. The 
reservoirs are in turn coupled to ground only through a 
capacitance C s t r - When only the four outer gates are 
energized, we say that the electrons in the 2DEG are 
confined in the "stripe" geometry. We observe signif- 
icant differences between the measured S-SET conduc- 
tance Gset versus applied gate voltage V g for the two 
different geometries, as will be discussed below. 

Regardless of the gate configuration used, we can ap- 
ply a single model of the environment to our results, as 
shown in the lower left in Fig. |]. The SET island is 
connected to its leads through junctions with resistance 



R\(2) and capacitance Cu2). We assume that the S-SET 
leads present an impedance Zg to the SET while the 
2DEG electrons have a total impedance Zm to ground 
which is coupled to the SET through the capacitance 
C*2D- Nearly the entire length of the SET leads is far 
from the Au gates, so that Zn is almost completely unaf- 
fected by the gate voltage V g . The electrons immediately 
beneath the SET are strongly affected by V g so that Z2D 
will in general be a function of V g , and may also depend 
on the configuration of gates used (i. e., on the pool ver- 
sus stripe geometry). Finally, the SET is also coupled 
to the Au gates by a capacitance C g . We neglect the 
possibility of a substantial impedance on the gate lines 
largely because C g m 20 aF is by far the smallest capaci- 
tance in the problem. Furthermore, any gate impedance 
would be substantially reduced since there are six gates 
whose impedances would combine in parallel. This gen- 
eral model (excluding the smaJl^ate capacitance C g ) has 
been investigated previously,EIrE3 but without consider- 
ing any particular form for the impedances Zg and .^d- 

In Table | below we give the relevant sample parame- 
ters for the two samples SI and S2 considered here. The 
parameters were determined from electrical measurement 
and simulations as discussed elsewhere£3 




FIG. 1: Schematic diagram of a typical sample, showing the 
Au gates, Ohmic 2DEG contacts, and SET leads and bonding 
pads. The upper right inset is an electron micrograph showing 
the S-SET island surrounded by the six Au gates. Application 
of a gate voltage V g causes QPCs to form at the locations 
shown. The S-SET is coupled to the 2DEG beneath it by 
a capacitance C20 as illustrated in the lower right inset. In 
the lower left we show a circuit diagram of our model of the 
S-SET environment, including the lead impedances Zt and 
the impedance Z2D associated with the electrons immediately 
beneath the S-SET, coupled to it through the capacitance 
Ctd- We also show the combined capacitance C g to the six 
Au gates. 
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TABLE I: Parameters for samples SI and S2. Capacitances 
are in aF and energies in fieV. 



sample 


Ci 


c 2 


c g 


C*2D 


E c 


E. h 


Ej 2 


SI 


181 


120 


20 


356 


118 


5.9 


3.5 


S2 


375 


260 


20 


382 


77 


27. 


16. 



2. Tunneling Rates 



In general for our samples the charging energy E c — 



e 2 /2Cs where 



+ C 2 + C2D + C satisfies 



E 



SS- A is the Josephson 



Ce — Ci 
j 3 < E c < fe B T, where £7^ — 
energy p£ junction j given by the Ambegaokar-Baratoff 
relatioro and A is the superconducting gap. Under these 
circumstances the S-SET island charge is well defined, so 
that charge states can..be used as the basis for calculating 
the tunneling rates.oE3 We will also be concerned with 
transport at sufficiently low bias voltages V and tem- 
peratures T, that we need only consider the tunneling 
of Cooper pairs, for which the-sequential tunneling rate 
through junction j is given bycj 



(1) 



which is valid for Ej j -C E c . Here 5fW = ff — f% is 
the change in free energy associated with the tunneling 
event, and the function P(E) describes the probability of 
the Cooper pair exchanging an energy E with the elec- 
tromagnetic environment during the tunneling process. 
Following the usual environmental theorycJ P{E) can be 
expressed in terms of the real part of the total impedance 
seen by the tunneling electron Re[Z t (w)] first through a 



J 



kernel K(t) 

/OO 7 
-00 ^ 



x{coth(- 



)[cos(wi) - 1] - zsin(wt)} (2) 



'2k B T 

and then through the Fourier transform 
1 f°° 

P{E) = — - / dt exp \K(t) + iEt/h] 
2nhJ_ 00 



(3) 



A calculation of the tunneling rate T must therefore be- 
gin with a clear understanding of the impedance Z t (io) 
presented to the tunneling electrons by the environment. 



3. Model of the Environment 

Given the circuit modeLshown in Fig. |l|, one can use 
standard network analysisEZI to calculate the impedance 
Z t {uj) seen by an electron tunneling through junction j 
in terms of the impedances aiid capacitances shown in 
Fig. 0. The result is given byl±J 



1 



where 



jC + Y 



c = 



Cji + C2D 



(4) 



(5) 



where j' = 2(1) for j = 1(2) and 



Y = 



(Cj, + c 2D )cyc 2 D + tLo{c r c 2D ) 2 (z e + z 2D ) 



{C r + C 2D )C / C , 2D [(Cj, + C 2D ) 2 + C),\Z t + C% B Z 2B + iuj(C y + C 2D )C r C 2 u(Zi + 2Z 2D )Z e 

I 



(6) 



ignoring terms of order C g /C 2 n- To proceed, we need 
accurate models of Zi and Z 2 d; we begin by considering 
Zt, 



4. Model of Z e 

Since our leads are fabricated above the 2DEG, which 
acts as a ground plane, it is appropriate to model them 
as transmission lines. lj The most general form for the 
impedance Z tl - of a lossy transmission line terminated in 
a load Zl is given by 



Z 



Zl + Zq tanli7£ 
Z{) + Zl tanh j£ 



(7) 



where Zq is the characteristic impedance of the line, 7 its 
complex propagation constant and i its length. At the 
relatively low frequencies (<• 10 11 Hz) considered here, 
it is reasonable to ignore the inductive reactance of the 
line and treat it as a simple RC line with a resistance 
and capacitance per unit length rg and eg. Doing so, we 
have that Zq = \J rg/iucg and 7 = y/iujcgrg. Looking out 
at the line from the sample, the line termination Zl is 
provided by the bias circuitry, which typically presents a 
low impedance < 50 fl. For simplicity we therefore take 
Zt, = in Eq. (0), and obtain the resulting approxima- 
tion 



re 
ituce 



tanh \J itorgcgi 2 



(8) 
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which we take as the basic form for the impedance of a 
finite RC line. This form has been considered previously 
in the context of incoherent tuimeling of Cooper pairs in 
individual Josephson junctionsE3 

While this is likely a fairly accurate description of the 
impedance of a section of our leads, when used in evalu- 
ating the kernel K(t) in (g) it leads to integrals which are 
analytically intractable. Fortunately a further simplifica- 
tion is possible. We are interested in the low energy part 
of P(E), which we expect from (|J) to be dominated by 
the long time behavior of K(t), which is in turn domi- 
nated by the low frequency part of the impedance Zrc ■ 
We therefore expand Eq. (ph around uj = to obtain 



i + (Lur eCe eyV6) 2 



(9) 



as a reasonable approximation to Zrc i n the interesting 
limit. | — | 

Another common treatmentc3 of the RC transmission 
line problem is to consider an infinite RC line, whose 
impedance is given by Z$ = ^Jrg/iujcg. Unlike the fi- 
nite RC line, for which the impedance Zrc approaches 
a constant rgl at uj = 0, the infinite RC line has a 1/\/uj 
singularity at u> = which dominates the long-time limit 
of K(t) and therefore P(E). The kernel K(t) for the 
infinite line, as well as P{E), can be calculated exactly 
in the T = limit. At non-zero temperatures, a feigh- 
temperature expansion must be performed instead.li3 



5. Model of Z 2D 

Having developed a model for Zg, we now consider a 
model for Z2D ■ The particular model will depend on the 
geometry we choose. For an unconfined 2DEG, the sim- 
plest choice is that Z 2 ^ is ohmic with an impedance re- 
lated to R sq of the 2DEG: Z 2D ~ R sq /3- When the elec- 
trons are confined in the pool geometry, they are coupled 
to the remaining 2DEG by two QPCs with conductance 
1/-Rqpc (assumed equal), which appear in parallel from 
the vantage point of the SET. There is likely some shunt 
capacitance Cqpc as well, but the associated roll-off fre- 
quency 1/-RqpcCqpc is typically large (~ 10 11 s _1 ) and 
we therefore neglect it. So for the pool geometry, we take 



Z2D — -Rqpc/2- 



(10) 



The stripe geometry is more complex. Here, in ad- 
dition to the QPC conductances, the electrons beneath 
the SET are coupled to two large electron reservoirs with 
resistance i? s t r located between the outermost Au gates. 
At their narrowest, the reservoirs are 0.6 /um wide, but 
broaden in five sections to a width of 500 /jm. Each sec- 
tion contributes roughly 2R sq , so that i? s t r ~ 10i? sq = 
200 These reservoirs are in turn coupled to ground 
capacitively through a capacitance C s t r , which we esti- 
mate from the size of the reservoirs to be on the or- 
der of 0.3 pF. Using Z^ = 2(i? Q p c + Z~ t l)~ x where 



Rstr + 1/ «wCstr we nn d 



Rc[Z : 



2D 



^QPC 



1 



2 T s 2 tr i? str /(-RQPC + Rett) 



,2^2 



where r str = C stT (R s t r + Rqpc)- For w 2 r 2 



(11) 
« 1, 



then, Re[Z2D] approaches i?Qpc/2, while for w 2 r 2 tr » 1, 
Re[Z 2 D] approaches \{Rs^ + R 



L QPCV 



The imaginary 

part of Z2D in the stripe geometry is nonnegligible only 
in the vicinity of uj ~ 1/T str , so for our purposes we ne- 
glect it. In general then, at low frequencies Re[^2D] is 
kept finite by the presence of the QPCs, and at higher 
frequencies is dominated by the smaller of -Rqpc/2 and 

-Rstr/2. 



6. Decomposition of Z t (u)) 

While the form for Z t (u>) given by Eqs. (Q)-@ is 
complete, it is generally too complex to make signifi- 
cant headway in calculating K(t). Fortunately, signif- 
icant simplification is possible. For typical values of 
re w lx 10 6 — 1 x 10 7 fl/m and typical line lengths I sa 0.5- 
1 mm, Zg 3> ^2D for small uj. In contrast, for sufficiently 
large uj, Zrc becomes quite small (< 10 fi) and the con- 
dition Z2D 3> Zg is usually satisfied. It then becomes 
possible to decompose Y in Eq. @ into a low u> part 
dominated by Zg and a high uj part dominated by Z 2 n. 

For small w, as long as Zg 3> Z 2 v, we can safely neglect 
the terms in Eq. involving Z 2 \y. Furthermore, for 
u) such that ui <C vam{l I CjZg, \jC 2 \yZg) we can ignore 
terms in Eqs. (^) and that depend explicitly on u>. 
Making these simplifications, we have that in the small 
u), large Zg limit 



Z t (u) « Y- 



C 2U ) 2 + Cl 



r v. - 



K\Ze 



(12) 



For large uj, we drop terms of order Zg/Z 2 D, and find 
that we can neglect the explicit frequency dependence 
in the denominator of rtq) for uj <§; cg/(Cj>) 2 ri ~ 1 x 
10 16 rad/s. In contrast, we cannot necessarily neglect 
the explicit frequency dependence in the numerator, and 
find 



Y 



— 

V C*2D / Z 2 D 



C\Cy 



(Cf+C 2B )C 



2D 



(13) 



We combine this with M) to find in this limit 



Z t (uj) 



1 



iuj(d + C 2 )<VC 2D + (Cs/Cao) 2 ^ 1 

TT^I- (14) 



iujCcff + k 2 Z 2D 

Combining this result with (|l2]), we obtain for the real 
part of Z t {uj) 



Re[Z t (uj)} = Kl Re[Zg(uj)} + - 

1 



K 2 Z; 



2D 



\ujC c ffK 2 Z 2D ] 



(15) 
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which we take as our basic model for the real part of the 
impedance seen by an S-SET fabricated above a 2DEG 
ground plane. We believe this model should be applicable 
not only to out, own system, but to that of the Berkeley 
group as well .113 




E (ueV) 

FIG. 2: (a) Re[Zt(u)] for three different environmental mod- 
els, all based on a transmission line with re = 2.9 x 10 6 id/m 
and c f = 1.0 x 10~ 8 F/m and £ = 6.94 x 1(T 4 m. For this 
graph we have used Z2D ~ 100 fl. Solid line: exact result 
based on the full forms for Z t (u>) and Ze- Dotted line: ap- 
proximate form based on the low frequency approximation to 
Zg and the decomposition of Re[Zt(w)]. Dashed line: infinite 
line result for the same values of re and ce- (b) Solid lines: 
P a (E) for the same transmission line parameters, for (top to 
bottom) T = 10, 20, 50, 100, 150, and 200 mK. Dashed line: 
P(E) for an infinite RC line at T = 0. 

To illustrate the degree of approximation associated 
with Eqs. (|) and ©, we show in Fig. |(a) Re[Z t {u))} 
for three separate models of the environment; all are 
based on an RC line with rg = 2.9 x 10 6 oj/m and 
eg = 1 x 10" 8 F/m. The solid line shows Re[Z t (w)] based 
on Eqs. and Eq. (|J), i. e., on the impedance of a 

finite RC line coupled to a ground plane with impedance 



Z2D, using the full form for Re[Z t (w)]. The dotted line 
is Re[Z t (ti;)] calculated using the low-frequency version 
of Zp given in (§) and using the decomposition ( |l~5| ) of 
Ke[Z t (oj)), while the dashed line is the impedance of an 
infinite RC line using the same values of m and eg (with 
no ground plane). We have not included a curve using 
the decomposition ([l5]) and the exact form for a finite 
RC line in (|8|) since it is virtually indistinguishable from 
the full Re[Z t (w)] shown. 

We note that the impedance of the infinite line rises 
above that of the more realistic forms for frequencies be- 
low a few MHz, so that in general it gives more weight to 
the low frequency modes and may be expected to give 
a more sharply peaked P{E). More importantly, the 
low- frequency approximation to Zt, while agreeing quite 
well below ~ 10 7 Hz with the exact finite line result, 
significantly underestimates it for intermediate frequen- 
cies below ~ 10 12 Hz. The approximation may there- 
fore be of limited use for larger bias voltages; for low 
biases, however, it is likely to be more accurate than a 
model based on an infinite transmission line, which over- 
estimates the impedance at low frequencies. Finally, at 
sufficiently high frequencies, the approximate and exact 
forms for Re[Z t (w)] converge. 

B. Calculation of K{t) and P(E) 

1. Calculation of K(t) 

Having produced a tractable form for Rc[Z t (w)], we 
can now proceed to a calculation of K(t) and P(E) 
through Eqs. (^) and (||). We begin by noting that when 
using the low frequency form for Z t in (^), both parts of 
Re[Z t (w)] have the same form, namely 

Re[Z t (uj)} = — — 2 2 (16) 

1 + UJ Z T Z 

for appropriate R and r. Calculations for K(t) an.d_ | P(t) 
for this form have been given in detail elsewhere, EJEil but 
emphasize a different range for R and result in different 
forms for P(E). Using 

1 R R Rt 2 uj 

10 1 + CJ Z T Z U) 1 + Id T 

and coth(tud/2k B T) = 1 + 2/[exp(?iw/fc B T) - 1] we find 
that 



K(t) = Rq 1 |jF c 





R~ 










Id 



2R 



exp(hto/kBT) 

-Tc. 



1 + ld 2 T 2 



dt — i sign(t) T t 
coth 



1 



R 



id 1 + UJ 2 T 2 
f hid \ 



\2k B Tj 



(18) 
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where T c [f{u))\ = 2 J '°° / (ui) cos Lot dco and T s {f(uu)] = 
2 J °° f(io) sinwi dui are the Fourier cosine and sine trans- 
forms of f{io) respectively, taken to be functions of 
\t\. We have ignored terms in K{t) independent of t; 
since for P(E) to satisfy the normalization condition 

JZo P ( E "> dE = 1 we must have K (°} = °i^ we wil1 later 
ensure normalization by adding an appropriate constant 

in any case. 

Of the four terms in curly braces in Eq. (p^jpEpr K (t), 
the first three can all be evaluated analytically.Ej An an=. 
alytic form for the entire kernel has also been found,c2l 
and analyzed for the overdamped case such that 1/r is 
large compared to the Josephson frequency ioj = nA/h. 
However, the range of R and t in which we are inter- 
ested was not investigated. Nevertheless, we have made 
some progress in certain limits. We note first that the 
fourth term in (18) depends on the temperature T only 
through the dimensionless combination T = Ti/ (2/cbTt), 
and write 



1 



■ coth 



2k B T 



= k T (t,T) 



(19) 



For zero temperature (T = oo), we have 



form may also be used for a sufficiently short and narrow 
section of transmission line. 

In the opposite limit, for which T < 1, we find that 
the approximate analytic result 



k T (t,T) w vre- |t|/T cotT 



(23) 



holds for times \t\/r > T. For typical temperature scales 
available in a dilution refrigerator (T = 20-400 mK), the 
logarithmic term is very small when— (p3|) is applicable. 
We can therefore neglect it, and findta that 



K s (t) = ~Lk B T\t\/h+^[cot(T) 



-isign(i)](e-l*l/ r - 1) 



(24) 



This result is typically most useful for dealing with fi- 
nite RC lines, for which r _1 is typically on the order of 
10 7 s- x . 

Finally, for 1 < T < 10, both the logarithmic term 
and fcr are of the same magnitude for the relevant time 
scales. The analytic forms for K(t) in ( p2j ) and ( |24] ) must 
then neglect some potentially important term. 



k-r(t, oo) = y/irGf 



2 1 



4r 2 





0,0, ■ 



(20) 



where G is a Meijer G function. In the long time limit, 
this result goes as —2(r/t) 2 . 

More generally, for T / 0, we find that it is impor- 
tant to consider the relative importance of the terms in 
Eq. (|l8|) . If we evaluate the integrals which we can treat 
analytically, we have 



2^1^^ 

r q { n 



In 2 



ln(l-e- 27rfeBT l*l/ s )+7 
\kr{t,T)\ (21) 



-l*l/n 



where 7 ~ 0.577216 is Euler's constant. In order to com- 
pare the relative size of the terms, we evaluate for numer- 
ically. We find that for T> 1, (either low temperature 
or small t), kq- decays slowly with time. For long times 
then the term going as e~^/ T is by far the smallest, and 
can be neglected. Of the remaining terms, in the long 
time limit the logarithmic term dominates over kq- and 
we write the kernel as 

Kl{t ) = _ 2i? + i n (i _ e-fcrfcaTltl/ft) 

Rq { Tl 



+ i~sign(i)-ln(7r/T)j (22) 

where the constant term ln(7r/T) will allow P(E) to be. 
approximately normalized. This is essentially the resultt2l 
of Wilhclm, et al, and is generally appropriate for deal- 
ing with the high-frequency part of Eq. (|l5|) due to the 
relatively small values of Z2D and C c ff. However, this 



2. Calculation of P(E) 

Having obtained analytic forms for K(t), a straightfor- 
ward application of Eq. (|^) allows one to calculate P(E). 
Letting g — Rq/R, we have from Eq. (|2|) for the large 
T limit 



Pi(E) 



(tt/T) 



2/ a 



2ir 2 k R T 



Re 



iE 



27rfc B T' 



1 



(25) 

where B(x,y) is, the beta function, in a gree ment with 
Wilhelm, et a^£3 For Z 2 d we have from ([15|) that R = 
K2Z2D and t = K2C c ff^2D- While (|2^) is only valid for 
g > 1, in terms of Z2D this condition becomes Z2D < 
Rq/k.2, so that the result remains valid for quite large 
Z2T> when K2 is small. |— ■ 
In the small T limit, we use (|24|) to obtairiiij 

P S (E) = -^-e' 3{T ^ ) Re\e-™ /9 1 2{'T,gy~ ll{T ' 9) 
x {r ( 7l (T, g j) - T ( 7l (T, g), 72 (T, g))} } (26) 

where T(x, y) is the incomplete gamma function, and 
71 = fT - if:, 72 = f (cotT - i) and 73 - |cotT. 
Typically, this will be applied to a finite RC line, for 
which R — K\rii and r = r/c?t 2 j 'v6- In some cases, for 
a short RC line, it may be more correct to use Pi(E), and 
substitute the appropriate forms for R and r in Eq. fl25| ) 
instead. 

We show P S (E) in Fig. |2|(b) for the same transmission 
line parameters as used in Fig. §(a). For comparison, 
we also show the T = form for an infinite transmis- 
sion line given by Pi n {(E) = eVo/2irE 3 e~ eV °/ 2E where 
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eV = (4r e n 1 /R Q )(e 2 /2c e ). Even at T = 10 mK, P s (£) is 
significantly broader than P in f(P). While it is possible to 
obtain an analytic form for P; n f(i?) by expanding Eq. (||) 
in the high temperature limit, the resulting expression is 
of limited use for E ^ 0, since it involves only even pop=. 
ers of E and therefore cannot satisfy detailed balance.Efl 
As a result, such an expression cannot be used to calcu- 
late I-V characteristics, for instance, whereas P S {E) in 
( p6| ) can. 

Ultimately, we are interested in calculating P tot (E) for 
the total impedance Re[Z t (cj)] seen by the tunneling elec- 
trons. If we were to calculate the total kernel K{t) for 
the decomposition in (p"5|), it would in general include 
all the terms in (|2l|). We were unable to find an ana- 
lytic form for Ptot(E) under those circumstances. How- 
ever, given the decomposition (|T5|), it is possible to write 
K(t) = K u (t) + K h{ (t) where Ku(t) and K h{ (t) corre- 
spond to the low- and high-frequency parts of Re[Z t (w)], 
with corresponding P\{{E) and P h f(E). The total P(E) is 
then given by the convolution P{E) = P\f{t) * P to t{E) = 
f-oo Pu(E - E')P M {E') dE', which can be performed nu- 
merically. 



C. Multisection Transmission Lines 

In our particular case, the sample leads do not have 
a single width. Instead, they broaden in sections from 
0.4 /zm (section 1) to 375 /im (section 4) as detailed in 
Table || below. As a result, we must generalize @ to 
allow for the possibility of multiple sections. In general, 
we use Eq. ([7) for a loaded transmission line, beginning 
closest to the SET with section 1. For this section, Zr, is 
taken to be the impedance of the second section, which 
is in turn terminated by the following sections. This cas- 
cading process is taken to end at our macroscopic con- 
tact pads, which are so broad as to provide very little 
impedance, and we therefore take Zr, = for the last 
section, so that its impedance is given by (H|). We also 
ignore a short (£ = 1 /im) section with w = 100 nm since 
it contributes only 50 f2 to Zi(0), and its associated P{E) 
is very sharply peaked around E = 0. 

If we were to use the exact form for Zg given by the 
above cascading procedure, it would be too complex to be 
of use. Fortunately, a simple approximation gives fairly 
accurate results. We take 



r c [z,H] = ]T- 



rgji 



(27) 



where rg i and cg i are the resistance and capacitance 
per unit length of section i, and £j is its length. We 
use the width Wi and length ii of each section along 
with the 2DEG sheet resistance R sq = 20 f2 and depth 
h = 50 nm to calculated rg i w R sq /{wi + 5.8ft.) and 
ct t ~ eeo(wi/h + 1.393), where e = 13 is the dielec- 
tric constant of GaAs. To find an approximate form for 
Re[Z t (u>)} we use the result (|^) for Zg in the decomposi- 
tion (Ely). For comparison, we plot both this approximate 



TABLE II: Transmission line parameters for the various sec- 
tions of the sample leads, with w and £ in fim, re in MQ/m, 
ce in nF/m, rgl in f2, and T calculated for T = 100 mK. 
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FIG. 3: Re[Zt(uj)] for four cascaded RC lines and a ground 
plane, using the parameters given in Table for (top to 
bottom) Z 2D = 6445, 3227, 1291, 258 and 100 fi. Solid lines: 
exact calculation of Re[/7t(u)]. Dashed lines: approximate 
version of Re[Z t (u)] as described in the text. 



result as well as the exact one obtained from the full form 
for Re[Zt(u)] in (f§)-(||) and repeated applications of (Q) 
versus frequency in Fig. ^ for different values of i^D- The 
agreement is very good, especially considering the num- 
ber of approximations required to develop a tractable 
approximate form for Re[^t(w)]. The approximate ver- 
sion tracks the exact result very well except between the 
various corner frequencies of the transmission line sec- 
tions, where its slope is generally too small. Agreement 
is better overall for larger values of Z2B > but even for the 
smallest values is still acceptable. 

We can calculate P\^ (E) for tunneling through junc- 
tion j for the four section transmission line by choosing 
either Pi(E) or P S (E) for a given section based on its 
value of T in Table O, and numerically convolving the 
four functions through 



P 



0) 



(E) = p(p(E)*P^(E) 



Pg(E)*P£>(E). (28) 



U), 



While somewhat time consuming, this procedure needs 
to be performed only once for a given temperature since 
the 2DEG beneath the transmission lines is not affected 
by the Au gates, and so the transmission line parameters 
do not change with V a . 

Finally, we then calculate the total P^l(E) for tun- 
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neling through junction j by convolving P^ (E) with 

P^2(E) calculated from Pi(E) for the appropriate value 
of Z 2 d- This procedure typically must be performed 
many times, but can be done relatively quickly. Re- 
sults for Ptot(E) for tunneling through junction 1 of S2 
are shown in Fig. |] for a series of different values of 
Z 2D . For Z 2B = 0, we take P${E) = P t {E), which 
is already relatively broad, with a width several micro- 
volts. In contrast, for small Z^d, -PodC-^) ^ s vei T sharply 
peaked around E — and approximates a delta func- 
tion, as can be seen in the insets (a) and (b) in Fig. ||. 

As a result, P^(E) is not strongly affected by P^(E) 
until Z2D ^ 200 ft. Finally, for sufficiently large Z 2 &, 

P 2 ^(E) begins to dominate and P to \ becomes very 
broad, indicating the high probability of inelastic tran- 
sitions. Overall, the trend is for the transmission line 
to dominate energy exchange for small Z 2 t>, while the 
2DEG dominates energy exchange for large Z 2 n. 



D. Calculation of I-V Curves 

To calculate the I-V characteristics for the S-SET, we 
use a master equation approachEj in which we assume 
that only two charge states, N and N + 1 where N is 
the number of Cooper pairs, are important. This should 
be a valid approach for temperatures and biases small 
compared to the charging energy E c of the S-SET. We 
begin by calculating the free energy change for changing 
the island charge from N to N + 1 (or vice versa) due to 
tunneling through junction j. We find 



80 



S f {l) 

°JN^N+l 



-Sf U) 

°JN+1- 

4E C {2N - 



■N 



+ 1) - (-l) j 2 aj eV (29) 



VgCg/e is the gate charge and no- 



where n g - 

(— l) J (Ci — C 2 )/2C% is the fraction of the bias voltage V 
appearing across junction j. We then use Eq. [j] to find 
the tunneling rates in terms of the sample parameters. 

The master equation can be solved exactly when only 
two charge states are considered. LJ Doing so, and using 
the detailed balance relation P{-E) = e- E/kBT P(E), 
we find 



I(V) 



ir 2 A 2 



sinh 



( eV 
\k B T 



32e 



Rk S cosh + S cosh ( 



2k B T 



(30) 



where Sf^ is the change in free energy for tunneling in 
the electrostatically favorable direction (N + 1 — * N for 
junction 1 and N — > N + 1 junction 2), fj — Rj/Rx, 



a/2 



the resistance quantum 



(=p5/0)), and R K = h/e 2 
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FIG. 4: Calculated P^(E) for S2, based on the transmission 
line parameters from Table |n] for a series of values of Z2D. 
Top to bottom: Z 2D = 0, 65, 129, 258, 430, 645, 860, 1291, 
1613, 2151, 3227, 4302 and 6445 fi. Insets: P^(E) for Z 2 u 
equal to (a) 65 (b) 430 (c) 1291 (d) 3227 and (e) 6445 Q. Note 
the scale change for Z 2 u = 65 and 430 fi. 



III. EXPERIMENTAL RESULTS 
A. Measurements 

We have performed electrical measurements on the two 
samples described in Table |in a dilution refrigerator at 
mixing chamber temperatures ranging from T = 20 to 
400 mK. Measurements were performed in a four-probe 
voltage biased configuration, with the bias applied sym- 
metrically with respect to ground, in a shielded room 
using battery powered amplifiers. High frequency noise 
was excluded with 7r-filters at room temperature and mi- 
crowave filters at the mixing chamber. 

Because the S-SET and 2DEG are electrically isolated 
from each other at dc, we were able to measure their 
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conductances Gset and G2D separately. In both cases 
the conductance was measured by applying an ac bias 
voltage at 11 Hz and measuring the resulting current us- 
ing standard lock-in techniques. The bias voltage used 
was 3 and 5 fiV respectively for the S-SET and 2DEG. 
We also performed measurements of dc I-V characteris- 
tics of the S-SET. For the 2DEG, the current contacts 
were positioned on opposite sides of the two QPCs, so 
that G2D measured the series combination of their con- 
ductances. Since the S-SET sees the QPCs in parallel, 
we take Z 2 n = 1/(4G 2 d) for the pool geometry. For 
the stripe geometry we cannot measure the stripe resis- 
tance i? s tr or capacitance G s t r directly, although we can 
estimate them from the sample design. We expect that 
in this geometry Z2B ~ -Rqpc/2 at low frequencies and 
Z2D ~ -Rstr/2 ~ 100 fl at high frequencies, as discussed 
above. While the I-V measurements of the S-SET were 
usually performed at a fixed gate voltage V g , the con- 
ductance measurements were typically performed versus 
V g for a variety of temperatures. In the pool geometry 
all six Au gates were tied together and V g swept from 
to —1 V, while in the stripe configuration only the four 
outermost gates were swept. 

Results of these measurements are shown for S2 in 
Fig. I for mixing chamber temperatures ranging from 
50 to 200 mK. As V g is made more negative, the S-SET 
conductance oscillates due to the effects of the Coulomb 
blockade; the oscillations are clearly visible in Fig. ||(a)- 
(e). We show results for the pool geometry for various 
temperatures in Fig. |5|(a)-(d), and results for the stripe 
geometry in Fig. |^(e). In both geometries, the envelope 
of the oscillations is relatively flat until V g ~ —0.3 V, 
at which point the amplitude of oscillations begins to in- 
crease. For the pool geometry, the envelope continues to 
rise until V g < —0.4 V, at which point it begins to fall 
again. For the lower temperatures (T < 100 mK), the 
drop with more negative V g is quite steep, whereas for 
the higher temperatures (100 mK < T < 200 mK) the 
drop is more gradual. For less negative V g , the envelope 
generally rises as T decreases, and tends to saturate be- 
low T w 100 mK. For more negative V g (beyond the 
maximum in the envelope at V g = —0.41 V) the envelope 
rises as T decreases until T w 100 mK, below which it 
decreases. 

Since G2D decreases as V g becomes more negative, 
the above behavior indicates that Gset varies non- 
monotonically as G2D decreases, first rising and then 
falling. A decrease in Gset as G2D decreases is expected. 
Physically, a larger G2D tends to damp phase fluctua- 
tions, promoting superconducting behavior and therefore 
a higher Gset- Alternatively, we can say that for low en- 
ergies, the probability P{E) of exchanging energy with 
the environment increases as G2D increases, as can be 
seen in Fig. [|. A higher G 2 d therefore implies a higher 
probability of elastic (or nearly elastic) transitions, and 
a higher Gset at low bias. In contrast, the decrease in 
Gset for V g > —0.4 V does not fit in with this general 
picture of energy exchange with the environment. While 
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FIG. 5: Gset vs. V g for S2 in the pool geometry for T — (a) 
50, (b) 100, (c) 150, and (d) 200 mK. In (e) we show Gset 
vs. V g in the stripe geometry for T — 25 mK. The inset to 
(d) shows an I-V characteristic for S2 at T = 25 mK with the 
2DEG left unconfined (V g = V). 



nonmonotonic behavior with G2D can be expectedjll it is 
generally associated with a crossover from the non-linear 
to linear portions of the I-V characteristic. We find that 
this is true in our simulations as well. 

In contrast, for our experiments Gset was always mea- 
sured in the linear part of the I-V characteristic. We 
show in the inset to Fig. ||(d) the I-V characteristic for 
S2 at T = 20 mK for V g = V, i. e., when the 2DEG 
is unconfined. The I-V characteristic is clearly linear to 
±8 /iV, so that our 3 fiV rms bias should be firmly in 
the linear regime. The non-monotonic behavior we ob- 
serve cannot then be associated with changes in P(E), 
and must arise from other physics. 

We can find a clue as to the source of this behavior by 
examining Gset vs. V g for the stripe geometry as shown 
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in Fig. Me). In this geometry the envelope also begins 
to rise at V g » —0.3 V. However, the rise is weaker, 
with the envelope increasing by only roughly half what is 
observed in the pool geometry. Furthermore, there is no 
decline in Gset as V g is made more negative. This last 
observation is in agreement with our model for Z 2 d in this 
geometry, which predicts that Z 2 d approaches i? s t r /2 at 
the higher frequencies which dominate our measurement 
of Gset- As a result, in this case P iot (E) never broadens 
for more negative V g as it does in the pool geometry, 
so that no decrease in Gset is observed. Once again, 
the environmental theory cannot explain the reduction 
in Gset for less negative V g . 

One possible explanation for this decrease is that 
there are changes in the offset charge of the S-SET is- 
land that occur on a time scale short in comparison to 
the time constant of the lock-in, but long in compar- 
ison to the time scales associated with environmental 
fluctuations. It is well known that charge fluctuations 
in the substrate give rise to 1// noise in SET-based 
electrometers. ocil Such charge noise typically has a mag- 
nitude of (Sq) 1 / 2 ~ 10- 4 -10~ 3 e/\/lFiat 10 Hz and a cut- 
off frequency (above which the intrinsic SET noise dom- 
inates) of about 100-1000 Hz. Let us assume that in our 
case the 1// noise is somewhat larger than is typical, say 
(Sq) 1/2 ~ 4x 10- 3 e/ at 10 Hz, due to the presence of 
the2DEG. If we write Sq = 1.6xl0" 4 e 2 /^ S /f, then 
the expected mean square charge variance^ between fre- 
quencies fi and f 2 is given by (<r 2 h ) = S ln(/ 2 //i). Tak- 
ing fx = 0.1 Hz and f 2 = 1000 Hz we find {(J 2 ch ) 1/2 ~ 
4 x 10 _2 e, so that a typical variance of a few hundredths 
of an electronic charge is not unreasonable. 

In a voltage biased configuration such as ours, these 
fluctuations would have the effect of averaging the mea- 
sured current over an ensemble of charge states centered 
around the gate charge n g . Similar effects hasffi been 
seen in measurements of other S-SET systems.c3 Since 
the S-SET current is sharply peaked around the charge 
degeneracy points, we expect that any such charge aver- 
aging would tend to reduce the measured peak current, 
and therefore the conductance Gset- It is therefore pos- 
sible that the reduction in Gset for less negative V g arises 
due to increased charge averaging as the electrons in the 
2DEG become less confined. We examine this possibility 
in more detail in the following section. 



B. Comparison with Theory 

To compare our measurements with theory, we must 
plot the measured Gset versus G 2 d- Having measured 
G2D versus V g , we can fit a smooth function to the mea- 
sured G2D and use it to convert V g to an approximate 
G2D- In Fig. [j|(b) below we show both the measured G2D 
and the fitted function versus V g in units of conductance 
quanta Go = e 2 /h. While the two are nearly indistin- 
guishable in the figure, this procedure in only useful for a 
limited voltage range. For V g > -0.31 V (G 2D Z 200G ) 
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FIG. 6: (a) Peak values of Gset versus G2Dfor T = 100 mK 
(A), 150 mK (O), and 200 mK (+). Calculated values G§ ET 
scaled to equal Gset at its maximum value at G£jf x at 200 mK 
are shown as the heavy solid lines, (b) On the left axis we plot 
measured G20 versus V g (solid line) , and the smooth function 
fitted to G20 vs. V g (heavy dashed line). On the right axis 
we plot values of a c h in the pool geometry for T — 100 mK 
(A), 150 mK (O), and 200 mK (+). We also show a ch for the 
stripe geometry (o). 



the measured values of G2D become unreliable. 

Having converted V g to G2D we must now plot the en- 
velope of Gset versus G2D; we do so by recording the 
positions of the peaks in Gset versus V g and using our 
fitting function to convert to G2D- The results of this pro- 
cess are shown for T — 100, 150 and 200 mK in Fig. f3](a). 
As can be ^Wrifl--Fig- ||j our data show no sign of 2e 
periodicity,0'E§K3c9 even at the lowest temperatures. A 
similar lack of 2e periodicity was observed by the Berke- 
ley group as well.tS At the present time, it is unknown 
whether this is an intrinsic feature of S-SET/2DEG sys- 
tems or is due to the inherent difficulties of excluding 
high-frequency-noise from all possible sources, including 
the substrateJid In any case, we presume that since our 
data are strictly e-periodic that probabilities of finding 
the S-SET island with cither an even or odd number of 
electrons are approximately equal. This will have little 
overall effect on our analysis, since at any given value 
of n g only one charge state contributes significantly to 
the current, the other being between charge degeneracy 
points and therefore in a low-current regime. Overall, we 
would then expect an additional reduction of a factor of 
roughly two in the measured current above and beyond 
that due to any other reason. 

In order to compare our results with theory, we calcu- 
late Ptot (E) according to the procedure described in Sec- 
tion II C above. We then use (|30|) above to obtain a cal- 
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culated conductance Gg Err . In general, the conductance 
calculated directly from (|3(]|) overestimates the measured 
conductance Gset by a significant factor (on the order of 
40-50). Such discrepan cies are not uncommon in small 
tunnel junction systems ,t§E§E3 but nevertheless compli- 
cate comparison with theory. The best we can achieve is 
to compare the relative change in Gset and G§ ET , but 
to do so requires that we pick some point of reference at 
which we will scale G SET to make it equal to Gset- 

Gset reaches its maximum value at some value of G^d 
which for all the temperatures considered here is approx- 
imately 6.5Go; we call this value G!, 1 D ? X , which may vary 
slightly with temperature. For G2D > G™^, Gset de- 
creases, behavior which cannot be accounted for through 
changes in the environment, making this region unsuit- 
able for choosing a reference point. For G2D < GJJJf*, we 
assume for simplicity that there is no charge averaging, 
and choose the point at which Gset reaches its maximum 
at T = 200 mK as our reference point. 

We plot G| ET versus G2D: choosing our reference point 
as discussed above, in Fig. g](a) as the heavy solid lines. 
All curves are scaled by the same factor (w 44) so that we 
can get some sense of the agreement in terms of tempera- 
ture dependence as well as dependence on Z2T>. Generally 
speaking, the agreement at T = 150 and 200 mK is quite 
good for G2D < G 2 1 g' x ; both calculated curves for Gg ET 
track the experimental values of Gset nearly exactly. We 
note that to obtain this level of agreement, we have used 
only one variable parameter, namely the scaling factor. 
All other parameters used in the theory are derived from 
experimentally measured quantities. Given the number 
of approximations involved, the level of agreement for 
these two temperatures is quite remarkable. 

For temperatures higher than T — 200 mK, the S-SET 
conductance begins to rise again. However, we expect 
the environmental theory discussed here to be applica- 
ble only for temperatures satisfying the Coulomb block- 
ade conditionka k B T < £ c i? Q /(27r 2 Re[Z t (0))], which for 
S2 corresponds to T w 240 mK for typical values of 
Re[Z t (0)]. Failure of the theory between T = 200 and 
250 mK is in good agreement with this condition. 

At T — 100 mK and below, the theory also disagrees 
with the experimental results. Specifically, the measured 
Gset does not rise as rapidly with decreasing T as pre- 
dicted by theory. Furthermore, the dependence of Gset 
on G2D changes from sublinear to superlinear, so that the 
S-SET conductance depends more strongly on the envi- 
ronmental impedance than theory predicts. This trend 
is accentuated at lower temperatures, as can be seen in 
Fig. ||(a) for T — 50 mK. In this case Gset is an even 
stronger function of G2D, and Gset is generally speak- 
ing slightly smaller that at T = 100 mK (not larger as 
would be expected from the theory). One possible ex- 
planation for a saturation of Gset would be that the 
electron temperature stops decreasing for some temper- 
ature below 100 mK. While it is likely that our electron 
temperature saturates (data at 20 mK differs only very 
slightly from the 50 mK data), such effects would not 



explain the change in dependence on G2D, or a decrease 
in Gset from 100 to 50 mK. 

There is also a lower temperature bound for appli- 
cability of the environmental theory, set by the condi- 

(i) (i) 

tion Pu 1 LxEj i <C 1 where Pmax is the maximum value of 

P^l(E). This condition must be satisfied for the per- 
turbative result for the tunneling rate in Eq. (Q) to be 
satisfied. In our case, we find that for tunneling through 
junction 1 (for which both Ej and P ma x are larger), 
PmL ~ 52, 41 and 34 meV" 1 , all at G 2D = 6.5G and 
at T = 100, 150 and 200 mK, respectively. In that case 
we find that PmaxE^ — 1.4, 1.1, and 0.92 for the same 
temperatures. In none of the cases is the condition for 
agreement with Eq. ([|) clearly satisfied, so that the agree- 
ment at 150 and 200 mK is perhaps better that might be 
expected. Still, if it were a failure of the perturbative ex- 
pansion which is leading to the disagreement at 100 mK, 
we would expect the theory to agree for sufficiently low 
G2D that P max drops to its largest value at 150 mK. In 
our case, this occurs at G2D ~ 4.0Go so that based on 
this argument we would expect theory and experiment 
to agree at 100 mK over much of the range shown, and 
only deviate for 4.0Go < G2D < 6.5Go- Clearly, this 
expectation does not hold for our data. 

We now turn our attention to the range G2D > G^g*, 
for which Gset decreases, in contradiction to the envi- 
ronmental theory. As discussed above, we consider the 
possibility that some form of charge averaging, which in- 
creases as G2D increases and the electrons in the 2DEG 
become more mobile, causes Gset to decline. To test the 
plausibility of this hypothesis, we calculate the average 
conductance (G SET ) given by 



(G 



SET/ 



w{n')G c SET (V,n')dn' 



(31) 



where G| ET (V, n) is the calculated SET conductance (in- 
cluding the scaling factor) calculated at bias V and gate 



charge n and w(n) = 



■ exp 



(,-i- 



iS assumed 



to be the probability of finding the SET in charge state 
n when the gate charge is n g . We vary <j c h to cause the 
(Gg ET ) to exactly match the measured conductance at a 
given V g and plot the results in Fig. ^(b) for T = 100, 
150, and 250 mK in the pool geometry, and a mixing 
chamber temperature of 20 mK for the stripe geometry 
(estimated electron temperature roughly 50-70 mK). By 
construction, <r c h = for G2D < G^f x in the pool ge- 
ometry. For the stripe geometry, a c h saturates at about 
2 x 10 _2 e for large negative V g , since in this geometry 
Gset never reaches as large a value as it does in the pool 
geometry. This result is physically reasonable, since the 
electrons in the stripe are never as confined as they are 
in the pool. As V g becomes less negative, a c h rises until 
it saturates (in all cases) at V g » —0.27 V, the voltage at 
which the 2DEG begins to deplete. The saturation value 
for <j c h is generally reasonable, being about 4 x 10 _2 e for 
the stripe and for the the pool at 100 mK, and 5 x 10~ 2 e 
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and 9 x l(T 2 e at 150 and 200 mK. While the cause 
of the increase in a c h at 200 mK is unclear, there does 
not appear to be significant further rise in <r c h for higher 
temperatures. 

To provide further support for this idea, we examine 
I-V characteristics for SI as shown in Fig. f!](a), which 
show the evolution of the I-V characteristics when the 
2DEG is increasingly confined. As the confinement is 
increased, the current initially rises (V g = —0.3 V) at 
all voltages, while the peak current remains at a fixed 
voltage. For Z 2 d = 1613 f2, the current has increased 
again, but the peak current has begun to move to higher 
bias. For Z 2 d = 2151 f2, the peak current has decreased 
and moved again to yet higher bias, while the current at 
higher voltages has generally begun to rise. Finally for 
Z 2 d = 6453 £1 the I-V characteristic has become quite 
broad and the peak has moved outward yet again. 

We can understand this evolution by examining the ef- 
fects of variations in Z 2 d and a c h on the I-V characteris- 
tics separately, as shown in Fig. ^(b) and (c) respectively. 
When Z 2 d alone is increased, the peak current drops 
and the voltage at which the peak occurs increases; at 
the same time, current at higher biases increases. This 
reflects broadening of P tQt {E) as Z 2 n is increased, and 
a higher probability of inelastic processes. In contrast, 
when <j c h alone is increased, the current decreases at all 
bias voltages, and the voltage at which the maximum 
current appears is more or less fixed. 

In order to obtain good agreement with experiment, we 
must include both variations in a c h and Z%d, as shown 
in Fig. @(d). Here, for V g = and -0.3 V, we take 
Z 2 d = and vary a c h, while for Z 2 \^ = 1613, 2151, and 
6453 fi we take a c h = 0. Overall the theory agrees with 
the experimental results quite well (apart from an over- 
all scaling factor), reproducing the initial rise in current 
with no shift in peak current position, followed by a re- 
duction in current and an outward shift in peak position. 
The agreement is poor only for Z 2 d = 6453 Q, for which 
the experimental current is too large in relative terms. 
Even here, however, the shape of the I-V curve is repro- 
duced nicely. Additional data (not shown) indicates that 
the S-SET current does usually drop sharply for large 
Z 2 £>, so that the amplitude of this particular I-V curve 
is probably anomalously high. 



C. Discussion 

Overall, the agreement between our experimental re- 
sults and calculations seems quite good, particularly for 
T = 150 and 200 mK. The good agreement of G§ ET 
and Gset for those temperatures, combined with the ac- 
curate predictions of our model for the evolution of the 
I-V characteristics gives us confidence that our model, 
despite its complexity, accurately describes our experi- 
mental system. In particular, it is clear that both the 
lead impedance and any impedance which is coupled di- 
rectly to the S-SET island must be included to give ac- 
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FIG. 7: (a) I-V characteristics for SI at an estimated electron 
temperature of 100 mK for an unconfined 2DEG (o), V g = 
-0.3 V (□), and Z m = 1613 (+), 2151 (o) and 6453 fi (x). 
(b) Calculated I-V characteristics (top to bottom at peak) 
for Z 2 n = 0, 129, 258, 430, 645, 860, 1613, 2151, 3227, 4302 
and 6453 Q,. Here a c h = for all curves, (c) Calculated I-V 
characteristics for (top to bottom) a c h = 0, 1, 2, 3, 4, 5, 6 
and 7 x 10" 2 e. Here Z 2D = for all curves, (d) Calculated 
I-V characteristics for an unconfined 2DEG (o), V g = —0.3 V 
(□), and Z 2 n = 1613 (+), 2151 (o) and 6453 Q (x). To fit 
the data at V g = and —0.3 V, we use a c h = 0.07 and 0.05e, 
respectively. For the remaining curves we take a c h — 0. 
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curate results. In our particular case, charge averaging 
appears to play an important role when the confinement 
of the 2DEG is reduced, either for less negative V g or 
when the stripe geometry is used. Overall, this improved 
understanding indicates that S-SET/2DEG systems can 
be used to test the accuracy of the standard environmen- 
tal theory in a way which was not previously possible. 

Generally speaking, the model of the environment pre- 
sented in ( [To] ) should be applicable to any system consist- 
ing of an S-SET and its leads fabricated above a 2DEG, 
including that of the Berkeley group. The model also 
gives a simple explanation for the lack of agreement be- 
tween the experimental results of the Berkeley group and 
the scaling theory of Wilhclm, et al. Since neither the 
leads nor the 2DEG can be ignored, P tot {E) must be 
a convolution of Pn{E) and P2d{E). In that case, our 
own calculations indicate that the power law behavior 
described in Wilhelm, et al., does not survive the con- 
volution. The resulting dependence of Gset on G*2D 
does however resemble power law behavior over a lim- 
ited range of G2D , with an exponent which can vary with 
temperature. Whether charge averaging effects are im- 
portant in the Berkeley system is unclear to us at this 
time. 

The most significant puzzle associated with our own 
work is the sudden disagreement between theory and ex- 
periment between T = 150 and 100 mK. While it is 
clear that the perturbative expression for the tunneling 
rates in Eq. [I] has a low temperature bound for applica- 
bility, it is unclear why theory and experiment disagree 
at 100 mK for all G2D < 6.5Go- This disagreement, 
as well as the continued reduction in current for tem- 
peratures below 100 mK is suggestive either of a limi- 
tation in the environmental theory, or that some other 
physics is beginning to dominate at the lowest tempera- 
tures. Note that while Z 2 d becomes relatively large, the 
impedance seen by the junctions never becomes much 
larger than about 1000 O at any frequency since both Ki 



and K2 are relatively small, so that we are still in the low 
impedance limit Re[Z t (w)] < Rq. Difficulties with the 
environmental theory in this impedance range could be 
of great importance to potential quantum computation 
applications. Finally, we do not believe that Coulomb 
blockade physics in the 2DEG is likely to be of impor- 
tance, since the QPC conductances are still quite high 
even for the smallest values of G2D; Coulomb blockade 
oscillations do not begin to appear in the 2DEG until 
both QPCs have a conductance below Go 

It is straightforward to propose experiments which 
could address these issues. The simplest would be to 
perform more measurements like those in Fig. ^| for sam- 
ples with smaller Ej. (Such measurements were unfortu- 
nately not performed for SI). One could also simplify the 
analysis significantly by constructing samples in which 
the 2DEG is selectively removed beneath the leads. In 
that case, the leads would no longer act as transmis- 
sion lines, and would present a small impedance to the 
tunneling electrons, so that in all probability only G2D 
would be of importance. We can also imagine reducing 
the coupling capacitance G2D so as to reduce the effective 
environmental impedance and keep the S-SET current 
relatively large even for large Z^d- This could be partic- 
ularly useful if the 2DEG is further confined so as to form 
a quantum dot, in which case the S-SET could be used 
to probe the dot impedance and energy level structure. 
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